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Preliminary  investigations  of  power  spectrum  estimates  for  signal  strength 
end  corresponding  saa-etate  time  fluctuations  are  reported  for  an  ovarwater  path 
between  two  drilling  platforms  in  the  Oulf  of  Mexico, 

Signal  strength  measurements  at  wavelengths  of  9*0,  5.3,  3*2  end  0/fb cms 
end  corresponding  ocean- wave  recordings  are  included  in  the  analysis. 

The  method  of  analysis  is  described  with  a discussion  of  important 
and  assumption  involved  in  the  techniques. 


I.  INTRODUCTION 


A recent  report  from  this  Laboratory  [1]  gars  the  results  of  the  reflection 
characteristics  of  centimeter  and  mill  lacier  radio  waves  for  a path  over  the  Oulf 
of  Mexico.  This  report  gave  these  characteristics  for  wavelengths  of  9.0#  5*3* 
3.2  and  0.86  cas  over  an  open-nee  path  between  two  drilling  platforms  off  Oread 
Island,  Louisiana,  (Figure  1), 

It  is  the  purpose  of  this  report  to  present  the  potter  spectre  estimates  of 
the  time  fluctuations  of  these  radio  signals  and  the  corresponding  ocean-wave 
recordings. 


XX*  DATA  INFORMATION 


The  details  of  the  radio  transmitters , receivers  and  ocean-wavs  recorder  ore 
described  in  Report  64  [1]  and  will  not  be  repeated  here.  The  transmitters  were 
located  on  the  seaward  drilling  platform  at  heights  of  either  13  ft.  or  33  ft  mal* 
The  receivers  were  on  the  shoreward  platform  at  heights  of  either  14  ft  or  53  ft 
msl.  The  platforms  were  5027  ft  apart  and  8 miles  off  shore.  Measurements  were 
made  with  the  antennas  both  vertically  and  ho ri tonally  polarised. 


A record  of  water-wave  height  vas  taken  with  a step  typo  wave  gage,  with 
steps  ev»ry  0.2  ft.  The  gage  was  located  near  a 14-in.  diameter  pile  on  the  sea- 
ward side  of  the  recoiver  platform. 
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III.  fcaATKEi  CC!  .71 
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I?.  DATA  SAMPLES 


Staples  of  the 
are  shown  on  figure 
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In  the  Appendix. 
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vys  of  the  tests  (6  and  7 of  August)  rain  showers  were  caasen. 

• southeasterly  flew  produced  relatively  low  amplitude  3 long 
X),  on  6 August , a squall  line  moved  into  the  area  from  the 
l gusts  up  to  33  ash  at  the  receiver  platform*  The  wind 
<st  through  ? August  with  an  average  speed  of  15  mph.  The 
id  with  wind  waves  averaging  13  to  24  in*  high  with  occasional 
i receiver  platform. 


. l data  taken  between  the  platforms  on  6 and  ? August 
It  be  pointed  out  that  these  samples  vers  not  taken 

: t the  water  data  are  merely  typical  of  the  sea  state  during 
presenting  a time  span  of  approximately  an  hour  and  a half* 
•.  he  data  for  each  run  at  a given  wavelength  were  made  with 
d the  sea  state  ds.ta  taken  at  the  same  time* 


DATA 


time  records  were  made  of  the  radio  aigaals  and  wave  heights 
-sites*  It  is  reeognicad  that  this  short  time  limits  the 
i.  Affects  of  this  limited  sample  of  data  are  discussed 


ial s were  received  and  recorded  with  a calibration  which 
v in  decibels.  The  proper  interpretation  of  a comparison 
■ , which  is  expressed  as  a logarithmic  ratio,  and  the  water 
sesd  as  a differ-'ffvce,  is  open  to  question. 

n!  the  wave  recorder  raises  the  question  as  to  whether  the 
re  of  the  *«a  state  over  the  path.  The  question  arises  as 
ghi  takes  at  a single  point  is  a good  representation  of 
uiity  of  the  pilings  of  the  drilling  platform  to  the  wave 
rtsd  to  affect  the  high-frequency  component 3 of  the  water 


i-sney  components  of  tr.e  power  spectra  are  further  subject 
5 Halted  frequency  response  of  the  SA  recorders  at  f ro- 
le per  second. 

~«suits  <euld  then  hinge  upon  the  following  assumptions 


nates  of  data  were  sufficient  to  insure  that  the  record 
rj  time  series. 
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TIME  VARIATIONS  IN  HORIZONTALLY  POLARIZED  RADIO  SIGNALS 
8ETWEEN  PLATFORMS  WITH  OCEAN  WAVE  RECORDINGS 
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2,  The  amplitudes  of  the  water  waves  for  frequencies  above  one  cycle 
per  second  are  negligible. 

3.  It  would  appoar  that  it  would  be  desirable  to  compare  the  sea  state 
with  the  signal  strength  in  7olts,  rather  than  in  decibels  so  that  each  of  the 
quantities  would  be  expressed  as  a difference  of  two  quantities.  The  labor  of 
ro plotting  the  original  data  to  a voltage  scale  and  the  errors  accumulated  in 
such  a process  prohibited  this  being  done.  Hence  it  is  assumed  that  the  frequency 
components  of  the  signal  strength  expressed  in  volts  would  be  the  same  as  when 
expressed  in  decibels.  This  has  been  checked  for  one  sample  of  data  involving 
the  3.2  cm  signal  (see  Figure  9),  where  this  was  found  to  be  a reasonable  assusption. 


VI.  ANmSIS  OF  THE  DATA 


The  power  spectra  of  the  data  were  made  by  first  obtaining  their  autocorre- 
lation function  with  the  correlation  computer  developed  at  this  Laboratory  [3]« 

This  function  was  analyzed  into  its  frequency  components  by  a method  presented 
by  Tukay  [2]  and  outlined  in  the  Appendix.  This  results  in  the  power  spectra  re- 
presenting the  average  power  in  a band  about  the  frequency  in  question. 

For  the  comparison  of  the  power  spectra  of  the  signal  strength  and  simul- 
taneous sea  state,  curves  are  plotted  on  a normalized  basis  (Figures  3*5,  8-11). 

In  order  to  show  the  relative  power  of  the  signal  strength  variations  for  the 
various  wavelengths,  the  power  spectra  of  these  variations  are  plotted  in  pro- 
portion to  their  rms  value  (Figures  6,  12).  A similar  ««t  of  curves  shows  the 
relative  power  for  the  sea  state  variations  (Figures  7,  13). 

i 

Since  it  is  of  interest  to  know  as  such  as  possible  about  the  amplitude  dis- 
tribution of  the  time  variations  when  dealing  with  correlation  and  power  spectrum 
analysis,  the  signal  strength  on  a decibel  basis  was  plotted  on  normal  distri- 
bution coordinates,  as  shown  in  Figure  14,  for  the  data  of  7 August.  This  was 
also  done  for  the  sea-stato  data  on  the  same  day  as  shown  in  Figure  15.  All  of 
the  power  spectrum  estimates  shown  were  for  vertically  polarized  radio  signals. 


VII.  RESULTS 


Figures  3.  4 and  5 show  the  normalized  power  spectra  at  the  radio  signals 
and  simultaneous  sea  state  for  wavelengths  of  0.86  cm,  3.2  cm  and  5.3  cm  taken 
an  6 August.  No  analysis  of  the  9-ca  signal  was  mads,  sines  the  signal  had  no 
appreciable  variation.  Figures  6 and  7 show  the  relative  power  spectra  of  these 
radio  signals  and  sea  states. 

Figures  3,  9,  10  and  11  show  the  normalized  power  spectra  at  0.66-,  3.2-, 

5*3-  and  9-cm  radio  signals  and  corresponding  sea  states  for  7 August.  Figures 
12  end  13  show  the  relative  spectra  of  this  data.  Figure  9 also  shows  the  results 
of  the  analysis  of  the  radio  signal  on  a voltage  basis  as  well  as  on  a decibel 
basis. 
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FIG  6-RELATIVE  POWER  SPECTRUM  ESTIMATES  (TUKEY  METHOD) 
FOR  53  CM,  3.2  CM,  AND  0.86  CM.  RADIO  SIGNALS 

AUGUST  6,  1952 


FIG.  7-RELATIVE  POWER  SPECTRUM  ESTIMATES  (TUKEY  METHOD) 
FOR  SEA  STATE  CORRESPONDING  TO  53  CM.,  32  CM 
AND  0.86  CM.  RADIO  SIGNALS  - AUGUST  6,1952. 
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FIG  8 - NORMALIZED  POWER  SPECTRUM  ESTIMATES  (TUKEY  METHOD) 
FOR  9.0  CM.  RADIO  SIGNALS  AND  SIMULTANEOUS 
SEA  STATE  - AUGUST  7.  1952.  TIME!  1125 
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FIG.  II  - NORMALIZED  POWER  SPECTRUM  ESTIMATES  (TUKEY  METHOD) 
FOR  086  CM  RADIO  SIGNALS  AND  SIMULTANEOUS 
SEA  STATE  -AUGUST  7,1952,  TIME'  1010 
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FIG  14  - ARITHMETIC  PROBABILITY  PLOT  FOR  RADIO  SIGNAL 
STRENGTH  IN  DB  FOR  WAVELENGTHS  OF  9.0, 

5 3,  32,  AND  0.86  CM.  ON  AUGUST  7,  1952 


FIG.  15 -ARITHMETIC  PROBABILITY  PLOT  OF  SEA  STATE  DATA 
CORRESPONDING  TO  RADIO  SIGNALS  AT  90., 

53,  3.2,  AND  086  CM.  ON  AUGUST  7,  1952 
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VIII.  CONCLUSIONS 


In  view  of  the  many  uncertainties  involved  and  because  cf  the  complexity 
of  the  phenomena  of  the  reflection  of  radio  signals  from  a rough  surface  any 
definite  conclusions  based  on  such  a 3hort  time  interval  and  so  few  samples  would 
be  doubtful.  It  is  felt  that  the  principal  value  of  the  information  is  that  it 
brings  out  points  which  scam  to  warrant  further  investigation  both  from  the  aKpari- 
mental  and  the  analytical  viewpoint. 

The  following  points  of  interest  may  be  noted: 

1.  The  correspondence  between  the  power  spectra  for  signal  strength 
and  sea  state  appears  to  be  better  for  the  longer  wavelengths. 

2.  The  power  spectra  for  the  0.86-cm  signal  strength  shows  little  simi- 
larity to  the  sea  state  and  contains  higher  frequencies  than  are  present  in  the 
sea  state. 


3.  The  rns  value  of  the  signal  strength  fluctuations  is  in  nearly 
every  case  larger  for  the  shorter  wavelengths. 

4.  Pcwer  spectra  for  the  sea  state  data  shows  considerable  variation 
between  runs  on  the  same  day  as  well  as  runs  taken  on  different  days.  This  would 
be  expected  since  the  weather  conditions  were  quite  variable  during  the  period  of 
the  tests. 


5,  The  signal  strength  distribution  is  approximately  log-normal,  but 
the  sea  state  data  shows  some  deviation  from  either  a normal  or  a log-normal  dis- 
tribution. 

As  a means  of  gauging  the  comparative  roughness  cf  the  sea  surface  to  dif- 
ferent signal  wavelengths  for  the  water  waves,  the  number  of  water  wave  peaks 
in  the  first  Frcsnell  zone  (see  Figure  16)  for  each  radio  signal  was  plotted 
against  the  frequency  of  the  water  waves  (Figure  17).  The  number  of  wave  crests 
was  plotted  from  the  equation  for  water-wave  length  for  wind  generated  [4]  in 
deep  wator* 


L . 5.12  T2 

where  L is  the  length  of  the  wave  in  feet  and  T is  the  period  of  the  wave.  Figure 
17  shows  that  for  shorter  radio  signal  wavelengths  there  are  fewer  wave  crests 
at  a particular  water-wave  frequency  than  for  longer- wav  e-length  radio  signals. 
If  the  fluctuations  of  the  radio  signals  v/ere  caused  by  the  horizonal  movement 
of  tho  wave3,  then  it  would  b8  expected  that  the  snorter  wavelength  radio  signals 
would  be  affected  by  the  water  waves  than  longer  wavelength  signals.  The  power 
spectra  of  the  signals  behaved  in  this  manner . 

To  investigate  the  offoct  of  the  vertical  movement  of  the  water  causing  the 
reflected  signal  to  pass  in  and  out  of  phase  with  the  direct  signals,  the  change 
in  water  height  required  for  such  a phase  change  was  determined  for  the  different 
signal  frequencies.  These  varied  from  0.5  ft  for  the  0.86-cm  signal  to  four  ft 
for  the  9.0  signal.  From  this  we  would  expect  the  shortsr-wave-length  radid 
signals  to  be  more  affected  by  the  vertical  movement  of  the  water  than  the  long- 
er-wave-length  signals. 
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APPENDIX 


Analysis  Techniques 


The  method  of  obtaining  the  power  spectirim  estimates  may  be  summarised  by 
the  following  steps: 

1.  Autocorrelations  were  made  on  the  original  data  by  a correlation 
computer  using  approximately  4000  sampling  points  for  40  time  lags.  In  each  case 
duplicate  autocorrelations  were  made  to  provide  a check  on  the  computer  operation. 
It  nay  be  stated,  however,  that  the  reliability  and  accuracy  of  the  computer  has 
been  well  established  in  previous  computations - 

2.  Prom,  a method  outlined  by  Tukey  [2],  the  “raw  estimates"  of  the 
power  spectrum  were  calculated  from  the  autocorrelation  by  performing  the  Fourier 
cosine  transformation. 
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where  An  is  the  autocorrelation  coefficient  for  nth  lag.  Dr.  Tukey  has  pointed 
out  that  these  “raw  estimates"  are  subject  to  serious  errors  because  of  the  math- 
ematical approximations  involved  and  the  propagation  of  substantial  statistical 
fluctuations  resulting  from  the  sampling  process  throughout  the  successive  steps 
in  the  utation. 

3,  Those  “raw  estimates"  may  then  be  smoothed  by  applying  the  correct 
factors  which  determine  the  interrelation  between  adjacent  values.  The  “smoothed 
power  spectrum  estimates''  are  thus  obtained  by  calculating 


Ufc  - 0.23  Si!_i  ♦ 0.54.  % ♦ 0.23  Sfc+1  (2) 

It  is  these  values  which  are  then  plotted  as  a function  of  frequency  in  Figures 
3 to  13.  This  results  in  the  total  power  between  frequencies  (h-l)X!  and  (h+l)J}, 
except  for  the  end  points. 

4.  Dr,  Tukey  has  further  shown  that  ’under  the  assumption  of  an  essen- 
tially Gaussian  distribution  for  the  original  data,  the  values  of  Un  are  distri- 
buted according  to  a chi-squared  distribution  with  f degrees  of  freedom. 
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where  K - total  number  of  sample  points 
m « total  number  of  lags. 
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Figure  Id  shows  the  behavior  of  a chi-squared  distribution  as  a function  of  the 
number  of  degrees  of  freedom.  Thus  for  50  degrees  of  freedom  the  estimate  will 
lie  between  0.74  and  1.45  of  the  true  value  9056  of  the  time.  The  accuracy  of  the 
estimate  increases  with  the  number  of  degrees  of  freedom.  This  increase  in  ac- 
curacy is  rather  gradual  above  50  degrees  of  freedom  but  drops  off  abruptly  as 
the  number  is  decreased. 

The  nunbor  of  degrees  of  freedom  for  the  curves  in  tnis  report  i3  approxi- 
mately 40.  Figure  19  illustrates  905b  limits  for  40  degrees  of  freedom  for  the 
5.3-cni  radio  data  on  6 August,  1952.  Included  on  this  curve  are  two  points  on 
the  90£  limit  curves  for  100  degrees  of  freedom.  This  would  correspond  to  about 
times  as  much  data  for  the  same  resolution  of  points  on  the  power  spectrum 
estimate. 

It  should  be  noted  that  for  a fixed  length  of  data  more  resolution  could 
be  obtained  only  at  the  expense  of  degrees  of  freedom  and,  hence, accuracy.  .In 
the  examples  in  this  report  approximately  40  degrees  of  freedom  were  chosen  to 
maintain  a given  probably  accuracy.  This,  of  course,  limited  the  bandwidth  or 
resolution  to  the  values  shown.  It  would  be  desirable  to  increase  both  the  re  — 
solution  and  degrees  of  freedom,  but  this  could  only  be  done  by  taking  longer 
samples  of  data.  It  may  be  stated  that  muen  longer  data  of  a similar  nature  is 
now  available  and  is  to  be  analyzed  in  the  near  future. 

It  is  also  very  important  that  the  number  of  samples  taken  from  the  auto- 
correlation to  obtain  the  "raw  estimates"  of  the  power  spectrum  be  sufficient  to 
detect  the  highest  frequency  components  in  the  original  data.  If  this  is  not 
done,  power  from  the  higher-frequency  components  will  be  aliased  into  (or  added 
into)  the  true  value  at  the  band  desired.  In  other  words,  power  actually  present 
in  the  record  at  higher  frequencies  will  appear  as  power  in  the  bandwidth  at  lover 
frequencies.  In  the  present  examples,  the  choice  of  sampling  interval  was  chosen 
to  detect  frequencies  up  to  about  10  cycles  per  second.  This  is  far  more  than 
necessary,  since  the  response  of  the  recording  instruments  is  limited  at  fre- 
quencies much  above  1 cycle  per  second.  This  choice,  however,  should  not  have 
limited  the  accuracy,  but,  instead,  simply  required  more  computation  than  may 
have  been  necessary. 

As  was  pointed  out  earlier,  this  approximation  of  probable  accuracy  i3  based 
upon  a Gaussian  distribution  of  the  original  data.  The  smallness  of  the  deviation 
from  a straight  line  variation  on  Figures  14  and  15  shows  that  this  was  approxi- 
mately true.  It  is  of  interest  to  note  that  the  signal  strength  expressed  in 
decibels  is  more  nearly  Gaussian  than  whan  expressed  in  volts.  It  was  found  that 
when  the  water  data  are  expressed  in  decibels,  the  result  was  no  more  Gaussianly 
distributed  than  when  expressed  as  the  original  difference  quantity. 
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